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Abstract 

We investigate the decay of an inhomogeneous Dl-brane wrapped on a with 
an electric field. The model that we consider consists of an array of tachyon kink 
and anti-kink with a constant electric flux. Beginning with an initially static con- 
figuration, we numerically evolve the tachyon field with some perturbations under 
a fixed boundary condition at diametrically opposite points on the circle 5"^. When 
the electric fiux is smaller than the critical value, the tachyon kink becomes unsta- 
ble; the tachyon field rolls down the potential, and the lower dimensional DO- and 
DO-brane become thin, which resembles the caustic formation known for this type 
of the system in the literature. For the supercritical values of the electric fiux, the 
tachyon kink remains stable. 
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1 Introduction 



Spatially inhomogeneous rolling tachyon on an unstable D-brane has been studied in 
boundary conformal field theory (BCFT) [I1I21I3], boundary string field theory (BSFT) [1], 
and Dirac-Born-Infeld (DBI)-type effective field theory [5l[6l[7l[8l[9l[l0l[IIl[l2l[T3l[Tl[T5]. 
In BCFT, this subject was considered in the presence of spacetime-dependent marginally 
deformed tachyon vertex operator at the worldsheet boundary. The evolution of the re- 
sulting energy-momentum tensor is qualitatively different from that of the homogeneous 
rolling tachyon [TOl [T7] . The energy density evolves into a localized delta-function ar- 
ray within a finite critical time. These delta-function singularities were interpreted as 
codimension-one D-branes [Il[2l[3]. 

In lower-energy effective field theory approaches for unstable D-branes, most of the 
actions involve only the first derivatives of fields on the brane. Since the actions are 
represented by truncating all the higher derivatives of fields, they are considered reliable 
only when the second and higher derivatives are small. Nonetheless, the effective actions 
have been shown to reproduce various nontrivial aspects of the unstable D-brane dynamics 
in special settings, for example, the spatially homogeneous rolling tachyon [181 CHI EDI EI] 
and lower dimensional D-branes as static worldvolume solitons [22l [23l [2H ESI EH EZl EH] ■ 

In an inhomogeneous time evolution of an unstable D-brane, however, its dynamics is 
not governed by a truncated effective action after a certain critical time. Felder et. al. 
showed that caustics develop in some regions during the evolution of the inhomogeneous 
tachyon field with a runaway potential in the DBI-type action [5]. They treated the 
tachyon as a coUisionless fiuid and found that the tachyon field becomes multi-valued in 
those regions within a finite critical time. They interpreted the caustics as a signal that 
the higher derivatives of the tachyon field blow up [5], [TT]. More examples for the pure 
tachyon evolution have been investigated in Refs. P, [H]. They observed that the slope 
of the tachyon field near the kink diverges within a finite time. Authors identified this 
phenomenon with the caustic formation observed in Ref. [5]. 

Interesting aspects of the low energy effective field theory with an electric fiux turned 
on, have been considered in the tachyon vacuum [29l EDI EH EDI EH [I5] . In this work, we 
shall examine the role of the electromagnetic fields on an unstable D-brane where those 
fields are varying in time and in space with a nonzero fiux. The electric fiux density here 
is identified with the string charge density on the D-brane [32l [33] . When the tachyon and 
the electromagnetic fields depend on time (or one spatial) coordinate only, the equations 
of motion imply that all the electromagnetic fields are constants [211 ES] . Consequently, 
the electromagnetic fields become simply parameters of the tachyon solution. However, 
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in an inhomogeneous time evolution process, the electromagnetic fields depend on both 
time and space in general [T5] . 

In this paper, we shall consider a simple case in the DBI-type action which contains 
a tachyon and a single electric field. The electric field is turned on along the inhomoge- 
neous direction which we restrict to be a circle S^. Assuming that the tachyon and the 
electric field depend on time and one spatial coordinate, we numerically solve the master 
equation of the tachyon field under a fixed boundary condition at diametrically opposite 
points, which is closely related to the formation of the DO-DO pair located at diametrically 
opposite points on the circle [SI [351 136] . 

We find that there exists a critical value of the string charge density under which 
the kink is unstable. The kink approaches a similar state which resembles the caustic 
formation observed in the pure tachyon case [8]. If the string charge density is increased 
above the critical value, we find that the tachyon kink becomes stable. 

In Section 2, we present the Dp-brane model with a tachyon and an electric field in 
the DBI-type action. In Section 3, we numerically solve the field equation and discuss the 
stability. In Section 4, we analyze the stability in a semi-analytic manner, and conclude 
in Section 5. 

2 D^^-brane Model with Tachyon and Electric Field 

We shall consider the dynamics of the tachyon field described by DBI-type action in the 
presence of gauge field interactions. With general gauge field interactions, the tachyon and 
the gauge fields are coupled in a very complicated manner. Therefore, we shall consider 
the simplest case in which the tachyon field depends on time and on one spatial coordinate 
only while the electric field is turned on along that inhomogeneous spatial direction. 

The DBI-type action with gauge field interactions on an unstable Dp-brane in flat 
space is given by [371 EH 139] 



where Tp is the unstable Dp-brane tension, V{T) is a runaway tachyon potential, and 
Ffj^i, is the field-strength tensor for the gauge field A^. The equations of motion for the 
tachyon and the gauge field read 




(2.1) 
(2.2) 




(2.3) 
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S.(3)^0, (2.4) 

where Cg^^^-j is the symmetric (anti-symmetric) part of the cofactor for the matrix (77 + 
OTdT + F)n^. We consider the tachyon and the electric field along the inhomogeneous 
direction living in the worldvolume of the Dp-brane, 

T = T(x°,xi), E = Foi{x\x'), (2.5) 

and, for simplicity, turn off all the other components of the electromagnetic field. For the 
unstable Dl-brane, the above field ansatze (12.51) are the general setting with no restriction. 
With these fields on the unstable Dp-brane, the quantity X and the cofactor matrix are 
simplified as 

X = -{l-f^ + T"^-E^), (2.6) 
(Cn = I ^/J^'^, I. (27) 

where T = OqT, T' = diT and all other components of C'^'^ are trivial, i.e., C'^'^ = 5^'^ . 
The solution to the gauge field equation (12.41) is then simply 

Ii = I3E = constant, (2.8) 

where 11 = dC/ d{doAi) = dC/dE is the conjugate momentum for the gauge field Ai in 
the Weyl gauge = 0. We also have a defined quantity 

(2.9) 

Using Eq. (12. 8p . we can express the electric field E with the tachyon field T and the 
constant electric flux IT as 

The tachyon field equation (12. 3p which completely governs the dynamics of the unstable 
D-brane, is now written as 

(1 + T")T - (1 - T')T" - 2TT'T' + " jp^^2y2 ^ = 0- (2-11) 

Setting n = in the above equation (12. lip , we obtain the field equation for the pure 
tachyon [5]. 
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The components of the energy-momentum tensor are 

r°o = (3{l + T'^) (2.12) 
= -pfr, (2.13) 
= -(3{l-f^). (2.14) 

Before we close this section, let us comment on what was studied for the pure tachyon 
case (n = 0) in a relation to a caustic formation. It is well-known that the tachyon po- 
tential for large T exponentially decays, V{T) ~ e~°''^ with an arbitrary positive constant 
a [ISj. Then for the case of pure tachyon field (11 = 0), the tachyon equation (12.111) is 
reduced to [5] 

(1 + T'^)T - (1 - f^)T" - 2fT'f' = a f 1 - + T'^) . (2.15) 



It was shown in the literature [HI El El [H] that Eq. (I2.15P develops a caustic at a finite 
time, which we shall reconsider in Sec. 4.3. 

On the other hand, with a nonvanishing electric flux 11, Eq. (12.111) is rewritten as 

(1 + T'^)T - (1 - f^)T" - 2fT'f' = -B{T){1 - + T'^), (2.16) 

where 

We shall discuss the role of the electric flux 11 for the 1/cosh-type tachyon potential in 
Section 4 by investigating the behaviors of B{T) for various II's. 



3 Evolution of Tachyon Kink Array on Dl-brane 

In this section we investigate the time evolution of an inhomogeneous tachyon kink array 
on an unstable Dl-brane in the presence of an electric field. The Dl-brane now contains 
a fundamental string with a string charge density 11. We shall also discuss the subsequent 
dynamical formation of the lower dimensional DO/DO-brane. Although we focus on the 
Dl-brane for simplicity, the results are generally applied to Dp-branes when only one 
spatial direction of them is inhomogeneous. We solve the time- dependent field equation 
(12. lip numerically with a well-known tachyon potential 



V{T) = 




(3.18) 
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For the constant electric field, E = Eq, with the above potential the static regular 



kink solution was obtained in Ref. [25]. For < -Eq < 1, the solution reads 



To{x) = A/2sinh 



— 1 sin 




(3.19) 



where the constant /3o = n/i?o is the value of (3 for the static solution. This solution 
represents a tachyon kink-anti-kink array on a Dl-brane, and reduces to a pure tachyon 
solution when 11 = |23j . 

In order to investigate the tachyon evolution, we evolve this static solution with a 
"minimal" perturbation by solving the time-dependent field equation ( 12. lip . Here, the 
minimal perturbation is implied by imposing the initial conditions, 

T{t = 0,x) = To(x), (3.20) 
f{t = 0,x) = 0. (3.21) 

In principle, if the the initial velocity of the tachyon field is zero, the system never evolves; 
it is easy to see from the field equation (12.111) that the static solution remains as the solu- 
tion in the course of evolution. However, in numerical simulations, the intrinsic numerical 
errors coming from the numerical scheme play a role of perturbation. Once the numeri- 
cal code is stable, the tiny numerical error induces a minimal perturbation to the static 
system. If the initial configuration is physically stable, the system will not evolve except 
exhibiting small fluctuations. On the other hand, if the system evolves from the initial 
configuration as time elapses, it means that the system is physically unstable. This is not 
from the numerical instability. Therefore, applying such a minimal perturbation is very 
useful for a system at present which is suspected unstable to form a caustic at later times. 
The boundary condition we impose is 

T{t,Xn) = 0, (3.22) 

where Xn is the location of the n-th node of the initial static configuration f l3.19p . This 
boundary condition implies that we fix the spatial coordinates of the kink/anti-kink cen- 
ters during evolution. By doing so, we may consider the time evolution of only one patch 
(a half period) of the kink array, and the rests are merely a duplication of that patch 
configuration. Under the boundary condition (3.22), the compactification radius is fixed 
by the periodicity of the static kink solution while it was not fixed in the Sen's work in 
Ref. [35]. 

The initial conditions, fl3.20p and fl3.2ip . imply the "point symmetry" about the kinks 
and the "reflection symmetry" about the line x = Xmax where the tachyon field has the 
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maximum value. The boundary condition fl3.22p that we impose preserves such symme- 
tries in the course of evolution. For the pure tachyon case, we can observe the formation 
of the DO-DO pair located at diametrically opposite points on the circle initiated by 
Sen, and with the given boundary condition the Z2 symmetry is preserved [SU [33, ES] • 
In the presence of electric flux, we impose the boundary condition in the same way as we 
do for the pure tachyon case, although the Z2 symmetry is now absent; we believe that 
it is an appropriate and the simplest set-up that we can consider in the first step for the 
stability check. 



3.1 Pure Tachyon Case (H = 0) 

Let us first consider the pure tachyon case (11 = 0) in which caustic phenomena have been 
previously observed in a bit different set-ups O [8]. In Ref. [5], the authors studied the 
caustic problem with an arbitrary initial configuration, and treated the tachyon field as 
collisionless fluid during the evolution. They found that there develop some regions where 
the tachyon field becomes multi- valued at a finite time. In Ref. [8], the authors studied 
the tachyon evolution in a very similar way to ours, but with different initial conditions 
and a potential V = exp(— aT^). Therefore, it is worthwhile to examine how the caustic 
phenomenon arises in our set-up if it exists. 

For n = 0, the static solution (13.191) is simplified as 



To{x) = y2sinh"^ 



(3.23) 



The proper range of is then given by 

0< Po< T^. (3.24) 

In the limit of /^g O5 ^^^^ slopes of the tachyon profile (13.231) at the kink and the anti-kink 
sites become infinite, and the region between the kink and the anti-kink goes to the closed 
string vacuum. However, the period of the static tachyon configuration remains at the 
same constant value 2\P1tx. Therefore, the solution becomes an array of step functions 
with an infinite height at each site of the kink and the anti-kink p2l [25] . For various 
values of /3o in the above range, we have performed numerical calculations. We observed 
that evolutions are more or less similar regardless of the value of /^o- We show numerical 
results in Fig. [1] which exhibit a typical evolution behavior. 

Beginning with the static configuration and a minimal perturbation, the tachyon field 
grows as it rolls down the potential. The velocity of the tachyon field T gradually grows 
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everywhere except the fixed nodes where lower dimensional DO/DO-branes are located. 
The initial static configuration is never physically stable. At late times, we can observe 
that the velocity near the kink grows faster than in the bulk, and the slope of the tachyon 
field gets very large. This behavior is very similar to what was observed in Ref. [5], so 
we identify it with the early stage of a caustic formation. When it evolves further, the 
numerical calculation diverges since the slope becomes too large as in Ref. [8]. 

Due to the caustic behavior at late times, it is not possible to observe the final decay 
products, such as lower dimensional D-branes and tachyon matter, in the inhomogeneous 
evolution of the unstable D-brane. Although the main purpose of present work is to 
investigate the stability of a given kink solution using numerical techniques, let us consider 
the decay products of our system in an analytic manner. Since we are considering the 
evolution of an unstable D-brane on a compactified circle, the total energy of the system 
is finite. The total energy of the initial kink (anti-kink) (13.231) is given by [23] , 



where the energy density T"*^ is defined in Eq. (I2.12p . The energy of the initial static 
kink (anti-kink) gives exactly the tension of DO (DO)-brane for the specific tachyon poten- 
tial (I3.18p . Note that the energy does not depend on /5o. In other words, the energy of the 
expected final decay product DO (DO) which corresponds to the kink (anti-kink) in the 
thin limit {(3q 0) is exactly the same with that of the kink (anti-kink) with an arbitrary 
(3o- If the story is so, no tachyon matter is produced in principle. For the other kink 
(anti-kink) solutions in more general runaway-type tachyon potentials [21], the analytic 
approach becomes more complicated. 

For our case, one can expect that the tachyon matter, if it should be produced in the 
very final stage, would come only from the minimal perturbation applied numerically, but 
it is a negligible amount. In the numerical evolution, however, such an ideal final stage 
can never be reached due to the caustic behavior as we described above, so the final decay 
products are beyond description. 



In this section, we consider the tachyon evolution with an electric field turned on. It is 
expected that the fundamental string which connects the lower dimensional DO/DO-branes 
modulates the tachyon fiow and decay. 

The electric fiux 11 = f]E remains constant during evolution. This is the solution to 
the gauge field equation as was discussed in Sec. 2. What is remaining is to solve the 




(3.25) 



3.2 Tachyon with Electric Flux {U ^ 0) 
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tachyon field equation, which we shall do numerically with applying 11 from a very small 
value. For 11 7^ 0, the static solution (13.191) is valid for 

f3^ - 7^ < < f^. (3.26) 

The value of /?o is not restricted, in principle, once 11 is ranged according to the above 
relation, or vice versa. 

3.2.1 Unstable Solution 

The pure tachyon solution (13.231) is a continuous limit (11 0) of the static solution (I3.19P 
with electric field. Therefore, turning on tiny amount of electric flux would not affect the 
evolution pattern so much, and the result for small 11 will be similar to that of the pure 
tachyon case. 

Since we are interested in the tachyon evolution in comparison with the pure tachyon 
case, first we fix the value of /^o and vary 11 in the range (I3.26p . When 11 is very small, 
we observe from numerical calculations that the tachyon evolves in a very similar way 
to the pure tachyon case. In Fig. [2], we show the numerical results for /5o = O.lTi and 
n = 10"3/5o (i.e., = 10^3). The tachyon field grows in time in the bulk, and around 
the kink/anti-kink at late times it exhibits a similar caustic behavior to what happened in 
the n = case. When the electric flux is small, the fundamental string cannot properly 
modulate the tachyon flow to the kink direction. This picture is common for all the values 
of /3o, and even for /3q > 1 which was a prohibited range for the pure tachyon case. 

As time elapses, the energy flows from the bulk region to the kink/anti-kink (See 
Fig. [3]). The energy density (12.121) in Hamiltonian representation is given by 



p = 7^ = y'(i + r'2)(n2 + n2. + y2)^ (3.27) 

where the conjugate momentum of the tachyon field T is 

dL 

Ht = ^ = (3.28) 

The plot of p = shows that the energy is lumped at the kink/anti-kink, and the 
corresponding lower dimensional DO/DO-branes become thin due to the nonvanishing 
momentum flow toward the kink location. As time elapses further, the evolution steps 
into the initial stage of the caustic formation in a very similar manner discussed in the 
pure tachyon case, beyond which numerical calculations cannot be performed. 

The pressure increases pretty homogeneously to zero from the below as shown in 
Fig. [3l The electric field which initially began from a constant value increases more or 



9 



less homogeneously in the bulk region, but it peaks up at the kink/anti-kink region at 
late times as shown in Fig. HI 

Similarly to the case of the pure tachyon in Section 3.1, we consider the decay products 
analytically. Substituting the static solution fl3.19l) into the energy density relation (12.121) . 
we obtain 



where ^ = -\/2/ \/ 1 — -Eg • Integrating over a half period centered at the kink, we obtain 
the initial total energy generated by the kink [25], 

^kink = / ' rfxT°° = .'^^^^^ + To, (3.30) 

where /3o = H/i^o and = v^ttTi. The terms on the right hand side of Eq. (13.301) denote 
the string charge over a half circle and the tension of the DO-brane respectively. The total 
energy of the kink profile depends on both 11 and /5o- The thickness of the kink can also 
be deduced from the slope of the static solution. 



To(x)~- JT,2-(/32-n2)x. (3.31) 
Po ^ 

For a given 11, once the tachyon is unstable and experiences changes in its thickness, the 
evolution would not be in such a way as to follow other static configurations which are 
described by different /5o's; the total energy -Ekink is not conserved. Being different from 
the pure tachyon case in this sense, it is hard to imagine what the final decay products 
would be. In particular, for our current system in which the electric field E evolves 
inhomogeneously near the kink, and which seems to accompany a caustic behavior at late 
times, it is impossible to discuss analytically the decay products at the final stage. 



3.2.2 Stable Solution 

As the value of 11 is increased, the electric flux on the fundamental string suppresses 
more efficiently the tachyon flow toward the kink as well as the tachyon decay in the 
bulk. Above some critical value of 11, therefore, we expect that the flux completely blocks 
the tachyon evolution and makes the tachyon stable. The existence of the electric flux 
(the T'^lP term) takes parts in the energy density (I3.27p . which presumably suppress the 
kinetic portion of the tachyon fleld for a given system. 
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From numerical calculations, we observe that the critical value ranges around Ilc{f3o) ~ 
O.iPo, or less. The numerical coefficient also depends on Pq. In Fig. [5l we show the 
numerical results for Pq = 0.1T[ and 11 = a/O.I/^q. The tachyon remains very close to the 
initial static solution with only small oscillations (these are visible in the zoom-in plot 
only). The velocity T oscillates about zero, which also indicates the stability of tachyon. 
In Fig. [6], we plotted the evolution of the maximum value of T (Tmax) for several values 

of n. 

For numerical calculations so far, we have adopted the initial condition fl3.2ip which 
implies only a minimal perturbation. Therefore, one may question what would happen if 
larger perturbations are applied to the stable solutions. The conclusion is that the stable 
solution remains still stable unless the perturbations are unreasonably large. We have 
performed numerical calculations with an initial condition, T{t = 0, x) = Ci x T(t = 0,x), 
by varying the constant Ci in a reasonable range (|ci| ^ 1). In Fig. [TJ we plotted the 
tachyon evolution for ci = 0.01 with the same values of the other parameters. The tachyon 
field remains still stable, although the oscillation amplitude gets larger than that of the 
minimal perturbation case. For larger ci's, the oscillation amplitudes are larger. 

According to our numerical results, the critical electric flux depends on the value of 
Po, HciPo) ~ O.iPo- For a given Pq, the stability of the kink is guaranteed if 11 > He- 
Conversely, for a given 11, Pq should range as 

n<Po< ion,. (3.32) 

The initial static kink with Pq beyond the above range, is subject to instability. The 
slope of the kink, ^JT^ — {Pq — U.'^)/Pq, also has a constrained range for the stability 
accordingly; the thickness of the lower dimensional D-brane does so, because it is evaluated 
roughly as the inverse of the slope. 



4 Semi- analytic Approach 

In the previous section, we observed from numerical calculations that the static tachyon 
configuration is stable when the electric flux is sufficiently large. For small or zero electric 
flux, however, the static configuration is no longer stable; the tachyon rolls down the 
potential and the caustic forms at the kink. In this section, we would like to capture 
the stabilization story more conceptually mainly by analyzing the field equation. The 
analysis will be a bit technical. 

The tachyon field equation (12. lip can be arranged as 

(1 + r^)f = (1 - t^) \t" + (1 + T'^) 5(T)1 + 2fT'f' + f^r^B{T), (4.33) 



11 



where 



^tanh(^), forn = 0, 



B{T)^{ r,^sech^(^).n.(^) (4.34) 
v/2[n2+r,2sech^(^)] ' ^ u. 

It is enough to consider only one patch between the kink {x = Xq = 0) and the anti- 
kink (x = Xi = -\/27r/3o/ \/13q— H^) as we did in numerical calculations. Let us call the 
location of the maximum of T as Xi/2 (= Xi/2). Equalizing the square-bracket part in 
Eq. f l4.33p to zero corresponds to the static equation. We adopted its solution as the initial 
configuration. Applying the minimal perturbation discussed in the previous section, we 
are now interested in how T behaves in order to check the stability. If its signature 
remains unchanged, T will grow/decrease only; the tachyon is unstable. In order to have 
a stably oscillating tachyon configuration, the signature of T should keep changing during 
the evolution. It is sufficient to look up this behavior only in some parts of the patch, so 
we shall do in the bulk region where the evolution is gentle. 



4.1 Pure Tachyon Case (E = 0) 

Suppose that the tachyon field begins to evolve initially with a tiny positive velocity 
T{t = 0,a;) > from the static configuration. The Noether momentum along the x- 
direction, 

Pi = -/5rr', (4.35) 

is always negative (positive) in the region < x < xi/2 (a^i/2 < x < xi). Therefore, the 
energy density fiows to the kink (anti-kink) accordingly, and there grows the slope of T as 
well. And the corresponding term, 2TT'T' + T'^T''^B{T), in Eq. ( 14.33^ is always positive. 

The remaining factor in Eq. f l4.33p in determining the change in T is now the square- 
bracket part 

r' + _ta.h(^)J. (4.36) 

Initially this is zero since we adopted the static solution as the initial condition. As time 
elapses, we observed from numerical calculations that the curvature T" remained more or 
less unchanged in negative in the bulk while T grows. Therefore, the above term fl4.36p 
remains positive since its second term is a growing function of T. 

As a whole, all the terms on the right hand side of Eq. fl4.33p remain positive, which 
keeps T positive; the tachyon is unstable to grow in the bulk region and to form caustic 
about the kink and the anti-kink. 
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The tachyon velocity T does not exceed one in the middle of the bulk region, which 
can be understood from the following. The slope T' and the curvature T" in this region 
is relatively small, so we can approximate the tachyon field equation to a homogeneous 
one 



-L(i-r^)t^h(^ 



(4.37) 



The solutions to this equation are given by [H 



v^sinh"^ Ac cosh (^-^ 
VSsinh^"^ AcCxp^-^ 
V2sinh^^ As sinh (^^j 



(4.38) 



where the constants Ac and Ag depend on the initial energy density and are less than 
while Ae is an arbitrary constant. Then |T| remains less than one for all of the solutions 
in Eq. 04.381) and approaches one at infinity. 



4.2 Tachyon with Electric Flux (U 0) 



If the electric flux is turned on, the only difference in the field equation fl4.33p is B{T) 
defined in Eq. fl4.34p . Unlike the pure tachyon case, B{T) is now a bump- like function of 
T of which maximum is at 

V(9n4 + 8n2Ti2)V2 _ 3n2" 



T = T 

± IT 



>/2sech" 



(4.39) 



We plotted B{T) for several values of 11 in Fig. [HI 
The square-bracket part 

T" + (1 + T'^ 



BiT) 



(4.40) 



which is initially set to zero can become negative if T is on the descending region of B(T). 
Then, it can possibly happen that the acceleration T becomes negative. If the decelerating 
state (T < 0) maintains, the velocity will change its signature and become negative, T < 0. 
Then T decreases and rolls back up the right hill of B{T). The square-bracket fl4.40p . T, 
and T evolve in an opposite way to that during descending. T will eventually change its 
direction in the end. This process keeps T oscillating about a certain configuration, and 
thus the tachyon becomes stable. The stabilization in this way is more probable for larger 
n since the maximum of B{T) is lower; a certain amount of electric flux is required in 
order to halt the tachyon rolling down the potential and the flow to the kink. This agrees 
with our numerical results obtained in the previous section. 
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4.3 Caustic Formation Revisited 



Before we close this section, let us consider the tachyon flow and the caustic formation 
in the kink region. The quantity —X in Eq. (12. 6p should be positive for a physical 
system. However, when a caustic formation is accompanied (11 < lie), it was observed 
from numerical calculations that the value of —X at the adjacent position to the kink falls 
to zero and then becomes negative. This behavior of —X occurs just before the numerical 
calculation diverges due to the large slope of T. This is very similar to the pure tachyon 
result in Ref. [8] where the authors used a different tachyon potential. 

Following the techniques in Refs. [5l [8], we can estimate the critical time at which 
the fluid element at the adjacent position reaches the kink. Once —X becomes zero at a 
position, assume that the tachyon evolution is constrained by 

_X = l_t2 + T'2-E2 = _^^(^l_T2 + r') =0, (4.41) 

where we used Eq. (I2.10p in the second step. As T^V'^/ {IP + T-^V'^) is nonzero, for 
n < lie, the system is governed by the equation 

1 - + T'^ = 0, (4.42) 

which is exactly the same one for the pure tachyon case. Then the results in Refs. [5], [8] 
are applied exactly in the same way as follows. 

This first order partial differential equation was solved analytically by the method of 
characteristics [5]. We introduce a parameter q which satisfies x{q,ti) = q, t{q,ti) = , 
and T{q,ti) = Ti{q) at the initial moment tj (we set = 0). Then the solutions to the 
equation (14.421) are 

Xe(g,t) = g-^S^t, (4.43) 
1 +T2 



T(q,t) = T.(q)+ =, (4.44) 

where Ti^g = dTi/dq and Xc{q,t) is the characteristic curve of q at t. After a finite 



time t = qyl + T^^/Ti^g, Xc will cross x = with a nonzero value of tachyon field, 
T{q,t) = Ti(q) + q/Ti^g. Finally the tachyon field becomes multi-valued there, and its 
slope diverges [8]. The caustic forms firstly at the minimum oi t = q\ h- + Tfq/Ti^g. 
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5 Conclusions 



We investigated the numerical evolution of an inhomogeneous Dl-brane with a string 
charge density. We considered a DBI-type action which contains a tachyon field T with a 
run-away potential = 1/ cosh(T/v^) and a gauge field Ai in a fiat (1+1) dimensional 
spacetime. The target space that we considered is a compactified circle with a fixed 
radius. The tachyon field presents an inhomogeneous kink and anti-kink array on the 
circle. The fundamental string which we considered induces an electric field E on the 
Dl-brane. The electric fiux 11 (the string charge density) is solved to be constant while 
the electric field is let to vary both in time and in space. 

There exists a static solution when the electric field is constant, E = Eq [231 125]. We 
focus our interest on the case of < -Eq < 1 for which the tachyon profiles represent a kink- 
anti-kink array. For the pure tachyon case {Eq = 0), it had been studied in Refs. [SI [H] 
that the tachyon field is unstable to form a caustic at the location of the kink, or the 
anti-kink; the tachyon field becomes double-valued at a finite time, and thus the lower 
dimensional D-brane becomes thin. 

For the initial configuration in evolving the tachyon field, we adopt the static solution 
of which configuration is determined by two parameters, the electric fiux 11 and the electric 
field Eq (or (3q = H/Eq). The compactification radius is then fixed to y/2/ a/1 — E^. The 
electric fiux 11 remains constant during the evolution as a solution to the gauge field 
equation, while E {(3) is allowed change. Applying a minimal perturbation on the initial 
configuration, we evolved the tachyon field. 

When the string charge density 11 is zero or sufficiently small, the tachyon is unstable 
regardless of the initial pressure of the tachyon kink, —(3q. T rolls down the potential V , 
and its slope gets large at the kink and the anti-kink. Although it is never possible to reach 
numerically the infinite-slope state, we expect that the tachyon configuration approaches 
the situation of caustic formation in a similar way discussed in Ref. [8]. According to 
the Sen's proposal in [1], the final state can be interpreted as a thin DO-DO pair located 
diametrically opposite points on the circle [Ml ESI l36] . 

When n is sufficiently large, the tachyon becomes stabilized. The initial static solu- 
tion remains stable even when larger perturbations are applied. The fundamental string 
connecting DO- and DO-branes modulates the tachyon fiow and provides the stability. The 
stabilization is achieved when the electric fiux is larger than some critical value 11 > He 
where He is a function of the initial /5o, and numerically it is Tic ~ 0.1/9o- The stable kink 
with n > lie is identified as the lower dimensional stable thick D-brane with thickness 
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It is not sufficient, however, to fully discuss the stability of the resulting thick DO-DO 
pair with our limited boundary and initial conditions. To clarify the properties of the 
system, one needs more investigations with various boundary and initial conditions. For 
example, if one imposes a deformed initial velocity of the tachyon and an appropriate 
boundary condition, one would be able to observe pair annihilation of kink and anti-kink 
which is related to the instability due to DO and DO charges. 

Another example will be to consider varying radius of in time, or to consider the 
evolution on Ri. In this case, the period of the tachyon kink which was initially set to 
2\f2T\ I a/1 — -Eq, may change in time. One possible way of evolution is to track the suc- 
cessive static solutions of different periods. The electric field then evolves homogeneously. 

Our numerical study in this work was the first step to investigating the instability of the 
tachyon-plus-electric flux system, with the conditions which preserve involved symmetries 
maximally. In the future, we hope our work is extended to other various directions. 
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Figure 1: Plot of (a) T{t,x) and (b) T{t,x) for the pure tachyon case with (3q = lO^^T^. 
From the bottom, the hues correspond to t = 0, 50, 100, 300, 310, 320, ...,370. (The hne 
style and color are in the same order in the rest of figures.) In (a), the tachyon field grows 
as it rolls down the run-away potential. Near the kink, the slope increases rapidly as the 
tachyon fluid flows to the kink. In (b), T increases gradually in the bulk region while it 
increases very rapidly near the kink. If it is evolved further, the slope gets very steep and 
a caustic is accompanied. 
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Figure 2: Plot of (a) T{t,x) and (b) T{t,x) for the tachyon-plus-electric field case 
with /3q = 10~'^T^ and 11^ = 10~^/5q. From the bottom, the lines correspond to 
t = 0, 50, 100, 300, 310, 320, 370, 374. The overall evolution is very similar to the 
pure tachyon case, but is lagged a little bit due to the nonvanishing electric flux. This 
case confronts a caustic in the end. 
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Figure 3: Plot of (a) energy density and (b) pressure T^^. As the tachyon evolves, 
the energy flows from the bulk to the kink and the kink becomes thin. (See the zoom-in 
plot of the kink region.) The pressure gradually increases to zero from the below. 
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Figure 4: Plot of electric field E. E grows more or less homogeneously in the bulk region, 
and peaks near the kink/anti-kink. 
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Figure 5: Plot of (a) T{t,x) and (b) T(t,x) for the tachyon-plus-electric field case with 
= 10-%^ and = lO-^Pl The lines correspond to t = 0, 30, 60, 600, but are 
distinguishable only in the zoom-in plot of the top region. The tachyon field is very 
stable and exhibits only very small oscillations. Although we have not plotted enough 
time slices, the top part of the tachyon configuration experiences about 20 oscillations till 
t = 600. In (b), T also oscillates about zero stably. (Depending on the values of i3q and 
n, T can be less wiggly.) 
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Figure 6: Plot of T^ax for Pl = IQ-^T-^^ and = 10-^ 10-^ lO'^, lO'^, IQ-^^l For 
subcritical values of 11, T^ax grows monotonically. For supercritical values of 11, T^ax 
becomes stable. 
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Figure 7: Plot of (a) T(t, x) and (b) T(t, x) for the case of a larger perturbation with 
an initial velocity T(0, x) = 0.01T(0, x). The other parameters are the same as in Fig. [51 
The tachyon field oscillates stably about the static solution as in the minimal perturbation 
case, but with a larger amphtude. 
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Figure 8: Profiles of B{T) for various 11. 
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